One-electron atomic-molecular ions containing 
Lithium in a strong magnetic field 



H Olivares Pilon^ *, D Baye^ 1^, A V Turbiner^ ^ and J C Lopez 

Vieyra^'^ 

1 Institute de Ciencias Nucleares, UNAM, Apartado Postal 70-543, Mexico DF, 
04510 Mexico 

^ Physique Quantique, CP 165/82, Universite Libre de Bruxelles, B 1050 Brussels, 
Belgium 

E-mail: *horop(§nucleares .unam.mx, dbaye@ulb.ac.be, 
^turbinerOnucleares .unam.mx, ^vieyraOnucleares . unam.mx 

Abstract. The one-electron Li-containing Coulomb systems of atomic type {li, e) and 
molecular type {li,li,e), {li,a,e) and {li,p,e) are studied in the presence of a strong 
magnetic field B < 10^ a.u. in the non-relativistic framework. They are considered 
at the Born-Oppenhcimcr approximation of zero order (infinitely massive centers) 
within the parallel configuration (molecular axis parallel to the magnetic field). The 
variational and Lagrange-mesh methods are employed in complement to each other. 
It is demonstrated that the molecular systems LiH'^'*', LiHe^'*' and Li2~'' can exist for 
sufficiently strong magnetic fields B > lO'' a.u. and that Lij'*' can even be stable at 
magnetic fields typical of magnetars. 
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1. Introduction 

The existence of very strong magnetic fields in the atmosphere of neutron stars 
B ~ 10^^ G (~ 10^ a. u.) and magnetars B ~ 10^^ G (~ 10^ a. u.) has posed different 
questions about the structure of matter under the influence of such extreme magnetic 
fields (see for a review e.g. p]). In particular, it raised the question about the possible 
atomic and molecular systems which can exist in a strong magnetic field: can molecular 
chains be formed? The interest in getting answers to these questions dramatically 
increased after the long-expected discoveries made by Chandra and XMM-Newton X- 
ray observatories of absorption features in the soft X-ray spectra of several magnetic 
stars. The first of that was the celebrated 2002 findings in the radiation coming from 
the isolated neutron star 1E1207.4-5209 of two wide absorption features at ~0.7 KeV 
and ~1.4 KeV [2]. One of possible explanations of these features might be the presence 
of the traditional and the exotic Hg^ [3] molecular ions in the atmosphere of this 
neutron star [4J. It is one of the reasons to assume that other traditional and exotic 
systems - atomic and/or molecular, containing one or more electrons - might be present 
in the atmosphere of the neutron stars and magnetars. It has recently been predicted 
that a variety of hydrogenic, helium and mixed hydrogen/helium molecular ions of one 
and two electrons can be formed when a sufficiently strong magnetic field is imposed 

[HIE]. 

It is already known that the simplest molecular ion exists for all magnetic 
fields, while two more one-electron exotic ions He2''' and HeH^'*' might begin to exist 
for magnetic fields B ^ 10^ a. u. and B ^ 10^ a. u., respectively (see [7|). Furthermore, 
for magnetic fields B > 10^ a.u. the exotic compound He2''' becomes the most bound 
one-electron system among the one-electron hydrogenic-helium ions. Then, a natural 
question arises: Can other one-electron two-center systems {Zi,Z2,e) exist in a strong 
magnetic field? A simple electrostatic analysis indicates that this is indeed possible for 
Zi,2 < 4 (see [5]). 

The goal of the present work is to study the existence of Lithium-containing ions 
11^+, LiH^+, LiHe^+ and Li^+ in the presence of a strong magnetic field. The relativistic 
corrections are assumed to be of a small importance for B < 10^ a.u. following the 
analysis by Duncan [8] that the longitudinal electronic motion is still deeply non- 
relativistic. The molecular axis is assumed to be aligned parallel to the magnetic field 
line since it is rather evident that this configuration is optimal if a magnetic field is 
sufficiently strong [5l [9]. The positively-charged centers are considered to be infinitely 
heavy (zero-order Born-Oppenheimer approximation). 

We use two methods. The variational method with physically relevant trial functions 
and the Lagrange-mesh method. The first one has been proven to be very efficient [5] 
with physically motivated trial functions. On the other hand, the Lagrange-mesh method 
[To] nowadays provides the most accurate results for the total energy of the ion. 
However, in practice, the implementation of the Lagrange-mesh method is made easier 

t 1 a.u. « 2.35 X 10^ G= 2.35 x 10^ T. 
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with an a priori knowledge of the equihbrium distance and the corresponding total 
energy. Our approach is the following: as a first step we use the variational method to 
obtain preliminary but already sufficiently accurate results for the equilibrium distance 
and the total energy and then the Lagrange-mesh method is used to check and possibly 
improve the variational results. 

Atomic units = e = h = 1 are used throughout, although the energy is given in 
Rydbergs. 



2. One-electron molecular ion at the Born-Oppenheimer approximation 

The Hamiltonian which describes two infinitely heavy centers of charges Zi and Z2 
situated along the z axis, and one electron placed in a uniform constant magnetic field 
directed along the z— axis, B = (0, 0, B) is given by 

^-^ .2 Z^ 2 Zo 2 Z-\ Zo , ^ , , ^ s , n / s 

H = -^ + -^-^ + {pA + Ap) + A^ , (1 

ri r2 R 

(for geometrical setting see figure [H the origin of the coordinate system is the middle 
point between the nuclei for all molecular ions considered here). The vector potential 
corresponding to a constant magnetic field B = (0, 0, B) is chosen in the symmetric 
gauge 

^ = |(-y,x, 0) . (2) 

Then the Hamiltonian takes the form 

.2 Z-i 2^2 2 Z-i Zn „ t. B^ n , . 

H = -A i ^ + + BL, + — p2 , 3 

ri r2 K 4 

where = + y"^ and is the z-component of the electron orbital momentum, which 
is a constant of motion. In practice , we consider the charges Zi = 3 and Z2 = 1, 2, 3 
for Lithium-containing ions. 

For the case of unequally charged centers the only integral of motion apart from 
the energy is the angular momentum projection on the magnetic field direction, and 
the eigenstates are labeled by a quantum number corresponding to the excitation, and 
a Greek letter a, vr, 5 corresponding to the magnetic quantum number m = 0, —1, —2, 
respectively. In the equilibrium configuration the eigenstates of the Hamiltonian ([3]) 
with equal charges {Zi = Z2) are characterized by two integrals of motion: and the 
spatial parity operator P [r ^ —r) with eigenvalues p = ±1. The parallel symmetric 
configuration is also characterized by the z-parity, Pz{z —>■ —z) with eigenvalues a = ±1. 
The magnetic quantum number m, spatial parity p and z-parity a are related by 

p = a(-l)'" . 

Thus, any eigenstate has two definite quantum numbers: m and p. Therefore the space 
of eigenstates is split into subspaces (sectors) each of them characterized by definite m 
and cr, or m and p. Notation involves the subscript g/u (gerade/ungerade) corresponding 
to positive/negative eigenvalues of the spatial parity operator P. 
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Figure 1. Geometrical setting for tlie system of one-electron and two charged centers 
Zi and Z2 placed in a magnetic field B in parallel configuration along the z-axis. R 
is the internuclear distance and ri, r2 are the distances between the electron and the 
charged centers Zi and Z2, respectively. 



3. Methods 

3.1. The variational method 

We first study the ground state of tlie Hamiltonian by means of tlie variational 
metfiod. We follow the general recipe for the choice of physically relevant trial functions 
presented in [H]. For the case of unequal charged centers we use a trial function of a 
similar form to the one used to study the HeH^'*' molecular ion in a strong magnetic 
field (see 0), 

trial = ^1 + ^2 + ^3 , (4) 



with 



i^i = Aie-'^^^'^e-^'^P' + Ase-^^'^e-'^^?''' , (5a) 

^3 = ^^e-^^^i-^^'^e"^*?"' + A5e-'''''-''^''e-^'^p' , (5c) 

where Ai,,.5, ai.,.8, Pi.. .5 are variational parameters. Considering the internuclear 
distance i? as a variational parameter, we have in total 19 variational parameters (a 
free normalization of the trial function (jlj) allows us to keep fixed one of the parameters 
A1...5). 

For the case of symmetric systems with equal charges Zi = Z2 = Z we use a ground 
state trial function which was successfully used to explore the and He2^ molecular 
ions in a strong magnetic field [5]. It has the same form (jl]) with permutationally- 
symmetric functions ( !5all5cl) . 

iji = Ai(e-"i'"i + e""^*^^) e-^i^^', (6a) 
^2 = A2 e-°2(-i+-2) e-P^Tp' , (66) 
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where Ai, ,3, ai.,.4, (3i...3 are variational parameters. Here we end up with ten variational 
parameters. The functions ?/'i^2 describe the coherent and incoherent interactions 
of the electron with the charged centers, respectively. Obviously, they are suitable 
modifications of the celebrated Heitler-London and Hund-MuUiken H^-trial functions 
by multiplication with the lowest Landau orbital. In turn, the function ip^ is a modified 
Guillemin-Zener function. It can be also considered as a non-linear superposition of the 
functions ■ipi^2- 



3.2. The Lagrange-mesh method 

With the aim to confirm and improve the variational results, we apply the Lagrange- 
mesh method p^fT2l[T3] which is an approximate variational calculation performed with 
the Lagrange basis and simplified by the use of the corresponding Gauss quadrature. 
So far, this method gives the most accurate total energy of H^-ion in a magnetic field. 
The implementation of this method is made easier with an a priori knowledge of the 
equilibrium configuration and the total energy of the system. Thus, the variational 
results obtained for all studied systems (/^,e), {li,li,e), {li,a,e) and {li,p,e) can serve 
as input information for the Lagrange-mesh method. 

In order to study a two-center molecular system of charges Zi and Z2 it is adequate 
to use the spheroidal coordinates {^,ri,ip) defined as 

(see figur^Il) with ^ G (0,oo), rj G (—1,1) and the azimuthal angle ip G (0,27r). Since 
the z-projection of the orbital angular momentum Lz commutes with the Hamiltonian 
([3]), it is replaced by its eigenvalue m (magnetic quantum number). Because we are 
interested by the ground state, the wave function should be nodeless, hence m = and 
it can be written as [10] 



*-=oW = p|yT7jV^o(e,r/). (8) 

Making a substitution of ([H]) into the Schrodinger equation for the Hamiltonian ([3]), we 
arrive at the differential equation for the ground state function ipoi^^v)^ 
4 



{T^ + T,) + V{^,V) 



with the effective potential 



v^o(e,^) = i^^o(e,^), (9) 



16 



^^e(e + 2)(l-r/2), (10) 



where 



J{i.v) = {i + lf-ri^ ■ (11) 
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Here 

We consider ipoi^yV) spanned in some finite basis 

^o(e,^) = X]5^Q,i^^,(e,^) , (13) 
i=i j=i 

where is the size in the ^-direction and Nrj is in the //-direction, respectively. 
Exphcitly, it is given by (see [TU] ) 



where 



/.(O = (-l)^(/.x.)V^^^e-^/^\ (15) 

4 - hxi 



9.W = (-l)""-'l/V^^. (16) 

2 rj — rjj 

here Jij = J{hxi,rjj) (see (ITT]) ), where Xi and r^^ are the ith zero of the Laguerre 
polynomial L]\f^{x) and the jth zero of the Legendre polynomial Pi\f^{T]), respectively. 
The dimensionless parameter h is introduced for convenience. Making its variation, it 
allows us to adapt the Lagrange mesh to the actual size of the molecular ion. Taking the 
wave function (fT3l) and using the Gaussian quadratures associated to each coordinate, 
the equation ([9]) gets the form of mesh equations (see [TO] ) 



y^}Tiji>j' + V{hxi, r]j) 6ij6i'j' - E 6ii'6jj']ci'jf = . (17) 
i'=i j'=i 

The kinetic energy matrix elements Tj^j/j/ are calculated in [TU]. The potential V{C,,ri) 
is evaluated at the zeros of the Laguerre (scaled by the dimensionless parameter h) 
and Legendre polynomials. Finally, the problem of solving the Schrodinger equation is 
reduced to searching eigenvalues of the matrix equation ([1] 



4. Results 



The results of the variational and mesh calculations for the m = ground state of the 
Lithium-containing one-electron Coulomb systems {li,e) as well as {li,li,e), {li,a,e) 
and {li,p,e) in magnetic fields B < 10^ a. u. in parallel configuration are presented in 
tables [D, [2], [3], and HI respectively. The results for the molecular systems, show that the 
corresponding potential energy curve starts to display a well pronounced minimum at 
finite internuclear distance at the threshold magnetic fields Bfh ~ 2 x 10^ , ~ 10^ and 
~ 10^ a.u. It provides a theoretical indication to the possible existence of the bound 
exotic diatomic molecular ions Li2^, LiHe^^, and LiH^"^. For each molecular system, 
the potential energy curve as a function of the internuclear distance R is characterized 
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by the presence of a potential barrier. At large internuclear distances the interaction of 
two charged centers becomes repulsive: the total energy curve approaches from above to 
the total energy of the atomic ion Li^^ (see Figs. El HI [7j). If the energy of Li^"*" is lower 
than the minimum on the potential curve, the system is metastable towards the decay 
to Li^^. All studied systems display two general properties of the Coulomb systems in a 
magnetic field: as the strength of the magnetic field increases, they become more bound 
(the binding energy grows) and more compact (the equilibrium distance decreases). 

In tables El [3l and H] we present the equilibrium internuclear distance Req, the 
corresponding total i^;™" and binding Ef, = B — ^g;™" energies, the dissociation energy 
Ediss = -E™*" — -E/"*^^ (corresponding to the decay into Li^"*"), as well as the position Rmax 
of the maximum E^"-^ of the barrier and its height i.e. the difference AE = ^™ax_^?Twn 
for the molecular ions Lig"^, LiHe^^ and LiH^^, respectively. 

For each molecular system we calculated the lowest vibrational energy E^^^ by using 
the harmonic oscillator approximation around the equilibrium position. It defines the 
zero point energy and is presented in tables El El and HI Finally, -Rcross indicates the 
value of the internuclear distance R beyond the maximum for which the total energy is 
equal to the energy value at the minimum. Req, Rmax, Rcross and AE give a qualitative 
description of the barrier. 

Variational calculations were performed in MacBook with an Intel Core 2 Duo CPU 
at 2.4 GHz. The minimization package MINUIT from CERN-LIB and an adaptive 
multidimensional numerical integration routine DOIFCF from NAG-LIB were used in 
the variational calculations. For the lowest-eigenvalue search in the mesh calculations 
the software code JADAMILU [14j was used. 

4.1. The atomic ion Li^'^ 

The hydrogenic ion Li^'*' exists for all magnetic fields and is stable. In table [1] we present 
the results obtained for the ground state total energy of Li^"*" in a magnetic field in the 
range 2.1 x 10^ < B < 10^ using both the variational method and the Lagrange-mesh 
method. In the variational approach, the trial function is taken as a linear combination 
of the products of the lowest Coulomb orbital by the lowest Landau orbital 



where .4, ai...4, /3i,..4 are variational parameters. This few parametric trial function 
provides highly accurate results for the total energy in comparison to the more accurate 
Lagrange mesh calculations (when converged) for the whole domain of magnetic fields 
studied (see table [1]). They are systematically higher than those obtained with 
the Lagrange mesh method. Their absolute accuracies increase from ~ 0.01 Ry at 
B = 2.1 X 10^ a.u. to ~ 0.04 Ry a.t B = 10^ a.u. In turn, the relative accuracy 
on the binding energies varies from 0.005% to 0.01%. A feature that arises here (and 
which is also found in other systems, see below) is the fact that when we increase the 
magnetic field, the number of points that are necessary in the Lagrange-mesh method 



4 




(18) 
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for a reliable (converged) calculation gradually increases as a magnetic field increases. 
Eventually, it limits the domain of application of the Lagrange-mesh method. For 
example, with available computer resources it was impossible to get a number of mesh 
points at -B = 5 X 10^ and higher magnetic fields in order to reach convergence. 

Table 1. Results for the total {Et) and binding (£'b) energies (in Ry) of the Isq ground 
state of the atomic ion Li^^ in a magnetic field using the variational method (upper 
lines) and the Lagrange-mesh method (lower lines) when available. 



B 


Et 


Eb 


2.1 X lO'^ 


20 825.331 
20 825.3232 


174.6695 
174.6768 


5 X 10* 


49 780.425 
49 780.4084 


219.5746 
219.5916 


10^ 


99 738.828 


261.1723 


99 738.7958 


261.2042 


5 X 10^ 


499620.917 
499620.883 


379.083 
379.117 


lO*' 


999 560.377 


439.623 


999 560.340 


439.660 


5 X 10^ 


4 999 395.71 
< 4 999 396.4 


604.29 
> 603.6 


10^ 


9 999 314.04 


685.96 



4-2. The molecular ion Lig 

The results of the variational and the Lagrange mesh calculations for the lag state of 
the one-electron system containing two Lithium nuclei, {li,li,e~) for magnetic fields 
B < 10'' a.u. in parallel configuration are presented in table [2l It is worth mentioning 
that for one-electron symmetric systems with two nuclei with charge Z at distance R, 
the total energy is related to the total energy of the by a relation similar to a scaling 
relation for the single nucleus case presented in p!5], 

E,{Z, B, R) = Z'E,{1, B/Z\ ZR) + ^^^^ ~ . (19) 

ri 

Though it is known that is stable at all studied magnetic fields < i? < 10'' a.u. 
(see [6]), the presence in the r.h.s. of the repulsive term 2Z{Z — 1) / R (which is absent in 
the single nucleus case) indicates that other one-electron molecular ions may not exist. 
For example, for Z = 2 the Heg"*" ion does not exist at i? < 10^ a.u. but it starts to exist 
for B > 10^ a.u. at first as a metastable state and then for B > 10^ a.u. as a stable 
one. For Z = 3 the Li2'^ ion does not exist for i? < 2.1 x 10^ a.u. (see the discussion 
below) . 
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Table 2. Molecular ion Lij^ in a magnetic field B: equilibrium intcrnuclear 
distance Req, total and binding Ef, = B — jrjmm g^gj-gigg jj^ (^j^g variational 

method (upper lines) and the Lagrange- mesh method (lower lines) ; dissociation energy 
Ediss = £'™™ — i?^'^^; Rmax of the maximum at the potential curve (see figure [5]); 
height of the barrier AE = E"^°-^ — E"^"^ and lowest vibrational energy Eq^''; Rcross 
is the distance for which Et = £'""". 



B 


o 

eq 


t 


Eb 




Edis 




J- XiYYi ft f 




Tpvib 





2.1 X 10^ 


911 

0.2118 


20 854.8373 


145 


1627 


9Q 

29 


uuo 
5141 


0.232 


0.0096 


n 1 1 

U. i i 

0.11 


n 947 
0.2118 


5 X 10^ 


0.139 
0.1394 


49 809.483 
49 809.3904 


190 
190 


517 
6096 


29 
28 


058 
9820 


0.252 
0.252 


2.90 
2.8843 


0.48 
0.48 


0.37 
0.3655 


10^ 


0.111 


99 764.908 


235 


092 


26 


080 


0.245 


7.99 


0.76 


0.45 


0.1107 


99 764.8113 


235 


1887 


26 


0155 


0.246 


7.98 


0.76 


0.4451 


5 X 10^ 


0.071 
0.0707 


499 628.776 
499 628.631 


371 
371 


224 
369 


7 
7 


859 
748 


0.225 
0.2246 


31.51 
31.6 


1.69 
1.69 


1.56 


10^ 


0.060 


999 554.225 


445 


775 


-6 


15 


0.215 


48.06 


2.33 




0.0597 


999 554.020 


445 


980 


-6 


32 






2.26 




10^ 


0.037 


9 999 225.689 


774 


311 


-88 


35 


0.186 


138.59 


5.98 





999610 




999550 I 1 1 1 1 

0.5 1 1.5 2 

R[a.u.] 



Figure 2. Total energy for the lug state of the Lij in parallel configuration at 
magnetic field B — 10^ a.u. as a function of the internuclear distance R. The Li^"*" 
energy (dotted line) and the lowest vibrational energy (dashed line) are displayed. At 
the minimum R^a ~ 0.06 a.u. 
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A trial function of form (jlj) is used in the variational calculations. For each value of 
the magnetic field shown in table [2l the upper line corresponds to the variational results 
while the lower line (when displayed) corresponds to the Lagrange-mesh results. For 
B = 10^ a.u., it was impossible to get a sufficiently large number of mesh points with 
our available computer resources in order to obtain a converged result. The converged 
Lagrange-mesh energies when available are better than the corresponding variational 
results. The absolute accuracies of the variational energies increase from about ~ 0.1 
Ry for B < 10^ a.u. to ~ 0.2 Ry at 5 = 10^ a.u. For the binding energy, the relative 
accuracy is ~ 0.07% at 5 = 2.1 x 10"^ a.u. while it is ~ 0.05% at B = 10^ a.u. 

For a magnetic field B ~ 2.1 x 10^ a.u. the potential energy curve starts to display 
a minimum for R = 0.21 a.u. indicating the existence of a metastable state of the 
molecular ion Lig^, unstable towards decay to Lig"*" — >■ Li^^ + li. Further calculations of 
the energy curve show that the potential barrier is not high enough to keep a vibrational 
level, i?Q** > AE. Hence, the system is unstable against vibrations. For magnetic fields 
-B > 5 X 10^ a.u., the potential energy well becomes sufficiently deep to keep at least 
one vibrational state. The height of the barrier increases when the magnetic field is 
increased and the system becomes more stable against vibrations, i.e., while both AE 
and E^^'' increase monotonously with the magnetic field, their ratio AE/Eq"^^ increases 
from ~ 6 for 5 = 5 X 10^ a.u. up to ~ 23 for B = 10^ a.u. At the same time, 
the dissociation energy E^iss = E"^"^^ — Ej'^'^ decreases and eventually, for magnetic 
fields B ~ 10^ a.u., the total energy becomes smaller than the total energy of the Li^"*" 
atomic ion. The system becomes stable towards the decay Lig^ — > Li^^ + li, a situation 
depicted in figure [2]for B = 10^ a.u. The qualitative behaviour of this molecular system 
is typical: when we increase the magnetic field strength, the molecular ion Lig"*" becomes 
more compact (the equilibrium distance decreases) and more bound (the binding energy 
increases) . 

The electronic distribution of the Lig"*" molecular ion in the lag state for 
B = 10^ a.u. is displayed in figure [3] for the internuclear distances Reg ~ 0.06, 
Rmax ~ 0.21, and R = 1 a.u. At equilibrium, the electronic distribution is characterized 
by two overlapping peaks centered at the positions of the charged centers. Eventually, 
at large internuclear distances the electronic distribution consists of two symmetric 
electronic distributions each of them corresponding to the atomic ion Li^"*" in a magnetic 
field. The schematic picture shown in figure[3]for the ground state electronic distribution 
of Li2^ is typical for all the magnetic fields where this molecular ion exists. 

4-. 3. The molecular ion LiHe^^ 

Table [3] contains the results of the variational and the Lagrange-mesh calculations for 
the la state of the one-electron system {li,a,e), for magnetic fields B < 10'' a.u. in 
parallel configuration. For this asymmetric system, a trial function of the form (jlj) 
is used in the variational calculations. For each value of the magnetic field shown in 
table [3l the upper line corresponds to the variational results while the lower line (when 
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Figure 3. Evolution of the electronic distribution J \ipo{x,y, z)\'^dy for the ground 
state lag of Li^^ in a magnetic field B = 10^ a.u. for internuclear distances: (a) 
Reg « 0.06, (b) R„,ax ~ 0.21 and (c) R = 1.0 a.u. 
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Table 3. Same as in table [2] for the molecular ion LiHe^^ in a magnetic field B. 



B 


Req 






Ediss 


Rmax 




Eq 


Rcross 


10^ 


0.110 


99 775.455 


224.545 


36.628 


0.144 


0.40 


0.56 


0.1663 


0.1101 


99 775.3783 


224.6217 


36.5825 


0.144 


0.3917 


0.55 


0.1661 


5 X 10^ 


0.067 
0.067 


499 654.950 
499 654.857 


345.050 
345.142 


34.033 
33.974 


0.136 
0.136 


9.229 
9.213 


1.72 
1.57 


0.2268 
0.2267 


10^ 


0.057 


999 589.934 


410.066 


29.557 


0.130 


16.630 


2.31 


0.2679 


0.0568 


999 589.810 


410.190 


29.470 


0.130 


16.610 


2.30 




5 X 10^ 


0.040 


4999 404.74 


595.26 


9.02 


0.117 


44.40 


4.26 


0.940 


10^ 


0.035 


9 999 309.05 


690.95 


-4.99 


0.111 


61.42 


5.35 





displayed) corresponds to the Lagrange-mesh results. For the binding energy, the relative 
improvement of the mesh calculations is ~ 0.03%, almost uniformly for magnetic fields 
B = 10^ — 10^ a.u. For i? = 5 x 10^ a.u. the Lagrange-mesh energy is worse than the 
variational one. 

For a magnetic field B ~ 10^ a.u., the potential energy curve starts to display a 
minimum for R ^ 0.11 a.u indicating the formation of a metastable state, unstable for 
the decay to LiHe^^ — * Li^^ + a. Also the potential barrier is small and it does not 
allow to keep a vibrational level Eq'^'' > AE. For magnetic fields -B > 5 x 10^ a.u., the 
potential energy well is sufficiently deep to keep more than one vibrational state. Also, 
the height of the barrier increases when the magnetic field is increased and the system 
becomes more stable against vibrations, i.e. the ratio AE/Eq^'' increases from ~ 5 for 
B = 5 X 10^ a.u. up to ~ 12 for B = 10'' a.u. At the same time, the dissociation energy 
Ediss = E^''"' — EY^"^ decreases and, eventually, for magnetic fields B ~ 10^ a.u., the 
total energy becomes smaller than the total energy of the atomic ion Li^^, indicating 
that the system can become stable towards the decay LiHe^^ Li^^ + a. However, 
even at 5 = 10^ a.u., the system is still probably metastable with respect to vibrations 
since E'q*^ > Ediss- The total energy Et as a function of the internuclear distance R for 
B = 10^ a.u., is shown in figure |H 

The electronic distribution of the molecular ion LiHe^^ in its Icr ground state 
for B = 10^ a.u. is displayed in figure [5] for three different internuclear distances: 
Req ~ 0.06, Rmax ~ 0.13 and Rcross ~ 0.27 a.u. The profile of the electronic 
distribution is asymmetric. At equilibrium, the electronic distribution for LiHe'^''' is 
characterized by one pronounced peak at the position of the Li nucleus and a shoulder 
due to the presence of the a particle. The electronic distribution evolves in such a way 
that eventually the shoulder disappears and the distribution becomes centered at the 
position of the Li nucleus and less asymmetrical. This ionic-bonding behaviour is typical 
for all the magnetic fields studied. 
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Figure 4. Total energy for the ground state Icr of the molecular ion LiHe placed 
in a magnetic field B = 10^ a.u. in parallel configuration. The Li^^ energy (dotted 
line) and the lowest vibrational energy (dashed line) are displayed. The equilibrium 
distance R^q « 0.057 a.u. 

In order to have a better picture of this evolution, figure [6] presents the ^-profile 
of the electronic distribution as a function of the internuclear distance for two values of 
the magnetic field B = 10^ a.u. and B = 5 x 10^ a.u. The electronic cloud becomes 
narrower in the 2;-direction when the magnetic field increases. 

4.4- The molecular ion LiH'^'*' 

In table H] we present the results of the variational and mesh calculations for the la 
state of system {li,p, e) for magnetic fields B < 10^ a.u. The trial functions (jll) are used 
for this asymmetric system. For a magnetic field B ~ 10^ a.u., the potential energy 
curve starts to display a minimum for R ^ 0.06 a.u. The potential barrier does not 
keep a vibrational level, Eq"^^ > AE. For i? > 5 x 10^ a.u. the potential energy well is 
sufficiently deep to keep at least one vibrational state. However, this system remains 
unstable towards the decay LiH^^ Li^^ + p for all magnetic fields studied. 

Only for the magnetic field B = 10^ a.u., we could obtain converged results with 
the mesh method. For the binding energy, the relative improvement was ~ 0.04%. For 
B = 5 X 10^ a.u., we also present the best Lagrange- mesh results obtained, although 
they are less accurate than the variational ones. 

A plot of the total energy Et as a function of the internuclear distance R for 
B = 5 X 10^ a.u. is shown in figure [71 The ground state electronic distribution for 
a magnetic field B = 10^ a.u. for three values of the internuclear distances: Req ~ 0.06, 
Rmax ~ 0.09 and Rcross ~ 0.12 a.u. is shown in figure [HI At equilibrium, the electronic 
distribution displays one peak which follows the position of the Li nucleus. Another view 
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Figure 5. Electronic distribution / IV'ola^,?/, z)\'^dy for the molecular ion LiHe^^ in a 
magnetic field B = 10^ a.u. as a function of the internuclear distance: (a) _Reg ~ 0.06, 
(b) Rmax ~ 0.13 and (c) Rcross ~ 0.27 a.u. 
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Figure 6. z-profile of the electronic distribution J J \iPq{x, y, z)\^dxdy oi the molecnlai 
ion LiHe"*^ for internuclear distances: Req, Rmax, R = 0.2 a.u. and magnetic fields: 
(a) B = 10^ a.u. and (b) S = 5 x 10^ a.u. 



Table 4. Same as table [2] for the molecular ion LiH^^ in a magnetic field B. 



B 






Eh 


Ediss 


Rmax 




rpvib 
Eq 


RcTOSS 


10^ 


0.0621 


999 587.579 


412.421 


27.20 


0.092 


0.95 


1.741 


0.1173 


0.0622 


999 587.42 


412.57 


27.08 


0.0928 


0.96 




0.1175 


5 X 10^ 


0.0415 
0.0415 


4999419.11 
< 4 999 419.3 


580.89 
> 580.7 


23.39 


0.089 


8.4 


4.403 


0.1668 


10^ 


0.0360 


9 999 333.65 


666.35 


19.61 


0.087 


14.2 


5.385 


0.2048 



is presented in figure [9] where the integrated electronic distribution along the 2;-axis is 
displayed for three different values of the internuclear distance and for i? = 10^ a.u. 
and B = 5 X 10^ a.u. In the latter case, the electronic distribution is almost symmetric 
around the position of the Li nucleus with a small asymmetry due to the presence of 
the proton. When we compare both fields, we can see how the electronic distribution is 
shrunk along the z-axis when the magnetic field is increased. 

5. Conclusions 

We presented a non-relativistic study of the ground state of the one-electron Lithium- 
containing ions Li2+, LiH^+, LiHe^+, Li^+ in the presence of a strong magnetic field 
B < 10'' a.u in parallel configuration within the Born-Oppenheimer approximation of 
zero order. Two methods were employed complementing each other: the variational 
method and the Lagrange-mesh method. 

The obtained results give clear indications that the exotic molecular ions Li2^, 
LiHe^^, LiH^"^ begin to exist as metastable states starting at the threshold magnetic 
fields Bth ~ 2 X 10^, ~ 10^ and ~ 10^ a.u., respectively. As the magnetic field 
increases the potential wells of all three systems become sufficiently deep to keep more 
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Figure 7. Total energy for the Itr-state of the LiK'^^-fon in parallel configuration at 
magnetic field B = 5 x 10^ a.u. as a function of the intcrnuclear distance R. The 
total energy of Li^^ (dotted line) and the lowest vibrational energy (dashed line) are 
displayed. The minimmn occurs at Req ~ 0.0415 a.u. 

than one (longitudinal) vibrational state. Eventually, the ions Lig"*", LiHe^^ become 
stable or almost stable at magnetic fields B ~ 10^ and B ~ 10'' a.u., respectively, 
when the ion LiH^'*' remains unstable towards decay LiH^'*' Li^"*" + p in the whole 
domain of magnetic fields considered B < 10^ a.u. The energies obtained with the 
variational method were improved with Lagrange-mesh calculations for magnetic fields 
B < 10^ a.u. for all three molecular ions studied (see tables [2], [3] and H]). For magnetic 
fields B > 10^ a.u., we were unable to reach converged results in the Lagrange-mesh 
calculations. 

The accuracy of the variational energy can be tested using the converged Lagrange- 
mesh results. For B ~ 10^ a.u., the relative accuracies on the binding energy are very 
close: ~ 0.05%, ~ 0.03%, ~ 0.04%, for Li^+, LiHe^+, LiH^+, respectively. 

All these molecular ions have a similar behaviour: when the magnetic field increases, 
each molecular ion becomes more compact (the equilibrium distance decreases) and 
more bound (the binding energy increases). For a given magnetic field where all three 
molecular ions display a minimum in the corresponding potential curves that support 
at least one vibrational level, we find a hierarchy of the binding energies, 

i.e., the most bound lithium-containing molecular system is Lig^. If we consider other 
ions with two centers and one electron [6j: He^+, ReR^^ and H+, in a magnetic field 
10^ < B < 10^ a.u., there is a hierarchy of ionization energies 




One-electron atomic-molecular ions containing Lithium in a strong magnetic field 17 



(a) 

10-^ 




Figure 8. Electronic distribution / \ipo{x,y, z)\'^dy of the ground state the molecular 
ion LiH^^ in a magnetic field B = 10^ a.u. for internuclear distances: (a) R « 0.06 , 
(b) R w 0.09 and (c) R w 0.12 a.u. 
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Figure 9. z-profile of the electronic distribution J J \ipQ{x,y, z)\'^dx dy of the ground 
state of the molecular ion LiH^"*" for internuclear distances: Req, Rmax, R = 0.2 a.u. 
and magnetic fields: (a) B = 10^ a.u. and (b) i? = 5 x 10^ a.u. 



The ionization energy grows with an increase of the total charge of the nuclei. In the 
case of the same charge, the higher binding energy corresponds to the presence of the 
nucleus with the higher charge. 
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